TOPICAL REVIEW 



Critical phenomena at perfect and non-perfect 
surfaces 

M Pleimling 

Institut fiir Theoretische Physik I, Universitat Erlangen-Niirnberg, D - 91058 
Erlangen, Germany 

E-mail: pleimStheoriel .physik.uni-erlangen.de 

Abstract. In the past perfect surfaces have been shown to yield a local critical 
behaviour that differs from the bulk critical behaviour. On the other hand surface 
defects, whether they are of natural origin or created artificially, are known to modify 
local quantities. It is therefore important to clarify whether these defects are relevant 
or irrelevant for the surface critical behaviour. 

The purpose of this review is two-fold. In the first part we summarise some of the 
important results on surface criticality at perfect surfaces. Special attention is thereby 
paid to new developments as for example the study of surface critical behaviour in 
systems with competing interactions or of surface critical dynamics. In the second 
part the effect of surface defects (presence of edges, steps, quenched randomness, lines 
of adatoms, regular geometric patterns) on local critical behaviour in semi-infinite 
systems and in thin films is discussed in detail. Whereas most of the defects commonly 
encountered are shown to be irrelevant, some notable exceptions are highlighted. It 
is shown furthermore that under certain circumstances non-universal local critical 
behaviour may be observed at surfaces. 
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1. Introduction 

Our current understanding of critical phenomena results from an intensive interplay 
between experimental studies of a large variety of physical systems, ground-breaking 
theoretical developments (including renormalization group methods, finite-size scaling 
theory and conformal invariance) and extensive numerical investigations of model 
systems. In many cases the systems under investigation are treated as bulk systems, 
thus neglecting the existence of surfaces which are unavoidable in real physical systems. 
Discarding surfaces in systems with short-range interactions is justifiable when studying 
bulk critical properties, as the contribution of the surface to extensive quantities is 
vanishing in the thermodynamic limit. However, a surface breaks the translation 
symmetry of a system and changes local quantities. Thirty years ago, it has been realized 
that this leads to surface critical behaviour which differs from the bulk critical behaviour. 
Since that time numerous theoretical and experimental studies have been undertaken 
in order to determine local critical quantities in systems with perfect surfaces. 

Real surfaces, however, are usually not perfectly smooth but display some degree 
of roughness due to the presence of surface defects as for example steps, islands, or 
vacancies. Impurities are also often encountered at crystalline surfaces and may be 
viewed as the source of some disorder at the surface. Furthermore, experimentalists 
nowadays create thin films which do not appear in nature, by growing artificial structures 
on the film surface. All these defects have some impact on magnetic surface quantities. 

The present work reviews the recent progress achieved in the study of critical 
phenomena in systems with boundaries. Besides semi-infinite systems and films with 
perfect surfaces, more complex geometries with wedges and corners as well as more 
realistic surfaces with defects are discussed. The review hereby focuses on the question 
whether the different types of geometries and/or of surface defects have an impact on the 
surface critical behaviour. It is therefore complementary to earlier reviews on surface 
criticality (HIS El El that exclusively considered fiat, perfect surfaces. 

The thermodynamics of a surface is completely described by the surface free energy 
per area Fg. Singularities occurring in Fg determine the phase diagram of semi-infinite 
systems. At some phase boundaries of this phase diagram surface and bulk free energies 
both exhibit singularities, whereas at other boundaries only Fg becomes singular. An 
example for the former case is the ordinary transition where bulk and surface ordering 
occur at the same temperature, whereas the latter case is encountered at the so- 
called surface transition where the surface layer alone orders, while the bulk remains 
disordered. This surface transition is encountered at temperatures higher than the bulk 
critical temperature, the critical fluctuations of the dimensional semi-infinite system 
then being essentially {d — l)-dimensional, corresponding to a phase transition in d — 1 
dimensions. 

In Section 2 the critical behaviour at perfect surfaces is discussed. At the bulk 
critical point different surface universality classes are obtained for every bulk universality 
class. These universality classes are discussed and their differences emphasized. 
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Furthermore, recent progress in our understanding of surface critical behaviour in 
systems with competing interactions is reviewed. Surface critical dynamics and the 
effect of symmetry breaking surface fields are also briefly discussed. Section 3 is 
devoted to more complex geometries with a wedge. Wedge-shaped models are very 
interesting systems where local critical exponents, which change continuously with the 
wedge opening angle, arise because of the particular geometrical properties of the wedge. 
However, at a given opening angle, surface critical behaviour at the ordinary transition 
is still universal and does not depend on microscopic details of the model as for example 
the lattice type or the strengths of the local interactions. This is completely different 
at the surface transition, where under certain circumstances non-universal local critical 
behaviour can be observed. At a fixed opening angle, edge critical exponents then 
not only depend continuously on the values of the local couplings but also reflect the 
existence of the disordered bulk. This intriguing behaviour results from the fact that 
at the surface transition the edge acts like a defect line in a two-dimensional critical 
system. Corner critical behaviour is discussed in this Section, too. In the last years a 
great deal of activity focused on the critical behaviour of wedges in presence of external 
surface fields. This rapidly developing field is also briefly summarized. 

Section 4 deals with the important issue of thin fllms and semi-inflnite systems 
with non-perfect surfaces. Surfaces arc very often naturally rough, due to the growth 
mechanism or because of erosion effects. Adatom islands, vacancy islands or steps 
are typical defects encountered at real surfaces. On the other hand, specific surface 
structures, as for example lines of adatoms or regular geometrical patterns, can be 
created on purpose by using advanced experimental methods. As these surface defects 
have an impact on local surface quantities, one has to ask the question whether they 
do change the surface critical behaviour. For semi-inflnite systems we must again 
distinguish between the ordinary transition and the surface transition. At the ordinary 
transition, common surface defects (presence of a step, surfaces with uncorrelated 
roughness, amorphous surface) are usually irrelevant for the STirfacc critical behaviour, 
but there are some notable exceptions. At the surface transition, defect structures 
like steps and additional lines of atoms located at the surface may yield non-universal 
local critical exponents, as observed in semi-infinite Ising models with additional surface 
structures. Interestingly, additional fines located at tfie surface of tfiin Ising films always 
lead to a non-universal local critical behaviour. Section 5 finally contains concluding 
remarks. 

2. Perfect surfaces 

The aim of this Section is to review the main results on surface criticality in systems 
with fiat, perfect surfaces. In Sections 2.1 and 2.2 we pay special attention to the 
phase diagrams of semi-infinite systems as well as to the sets of critical exponents 
characterizing the different surface phase transitions. As the amount of papers published 
on this topic is rather impressive, only selected results are presented. The interested 
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reader is referred to the reviews of Binder p , Abraham [2] , and Diehl j^l H] for further 
readings. The main experimental techniques available for the investigation of surface 
critical behaviour are reviewed in the book of Dosch [H]. Section 2.3 is devoted to critical 
phenomena in thin films, whereas in Section 2.4 the recent studies of surface criticality 
in systems with competing interactions are discussed. Brief accounts of works on surface 
critical dynamics and on surface critical behaviour in presence of external fields close 
this Section. 



2.1. Surface quantities and phase diagrams 

When comparing systems with and without surfaces, one remarks that extensive 
quantities are altered by the presence of surfaces. Thus, the free energy is given by 

F = F,V + F.,S + ... (1) 

with the bulk free energy per volume Ft, and the surface free energy per surface area Fg. 
V and S are the volume and the surface of the system, respectively. Further terms in 
equation ((H) may result from the presence of edges, as discussed in Section 3. Similar 
corrections also show up in other extensive quantities. 

The surface free energy Fg completely describes the thermodynamics of the surface. 
Similar to Fb, singularities occur in Fg. The values of the variables of Fg, where these 
singularities appear, determine the boundaries separating the different phases of the 
surface phase diagram. As we shall see later in the course of this Section, at some phase 
boundaries both F^ and Fg exhibit singularities whereas at other ones only Fg becomes 
singular. 

Local quantities as the local magnetisation density (often called local magnetisation 
for brevity) or the local energy density are also modified by the presence of a surface. 
Examples of order parameter profiles obtained by Monte Carlo simulations of three- 
dimensional Ising films with L = 80 layers at temperatures below the bulk critical 
temperature ksTd Jb = 4.5115 are shown in Figure^ [Hj. Here, all the couplings involved 
have been chosen to have equal strength Jb and no magnetic fields have been retained. 
The surface magnetisation mi = m(l) = m{L), related to Fg by the equation 

OF 

- -Ml' 

is smaller than the bulk magnetisation due to the reduced coordination number at the 
surface. The field Hi acts exclusively on the surface layer. The magnetisation increases 
from its surface value to the bulk value at distances exceeding the bulk correlation length 
^b- In case the thickness of the film is not large compared to E,b (as for example at the 
temperature ksT/Jb = 4.49 in Figure the bulk magnetisation is never reached. For 
surface couplings exceeding the bulk couplings by a large amount, magnetisation profiles 
decreasing monotonically towards the center of the system may also be observed. 




Figure 1. Order parameter profiles m{z) of a tliree-dimensional Ising film with 80 
layers at three different temperatures ksT/Ji,. The dashed lines denote the bulk values. 
After 0. 



For Ising models in film geometry we have the general Hamiltonian 

H = -Jb^SiSj - Js SiSj -i^^Si - Hi ^i- (3) 

(ij) (ij) i surface 

surface 

where the spins Si can take on the values ±1. The first sum runs over bonds connecting 
neighbouring spins where at least one of the spins is a bulk spin, whereas the second 
sum runs over all surface links. Bulk couplings, Jb, and surface couplings, Jg, are 
ferromagnetic. The bulk field H acts on all the spins, the surface field Hi only on spins 
located at the surface, i.e. a surface spin sees the field H + Hi. Sometimes a different 
coupling constant is used for the bonds connecting surface and bulk spins jHI- In general, 
the perturbations due to the presence of a surface are supposed to be of short range. 

Analytical results on surface critical behaviour are commonly obtained in the 
framework of continuum field theory, see [3 E] for comprehensive reviews. The standard 
0^ Ginzburg-Landau model, appropriate for the universahty class of the Ising model, is 
thereby augmented in order to include contributions coming from the surface: 

fs:=Fs/kBT=^co(l)^-hi<f) (4) 

with hi := Hi/ksT, whereas Cq is related to the surface enhancement of the spin-spin 
coupling constant in the corresponding lattice model. The resulting Ginzburg-Landau 
free energy density is readily generalized to cases where the order parameter exhibits a 
different symmetry. 

Layer- dependent quantities are very useful when investigating systems with 
surfaces. Examples are the magnetisation per layer m{z) and the susceptibility per 
layer x{^) where z labels the layers parallel to the surface. The surface magnetisation 
is nil = = 1) and the local susceptibility at the surface, i.e. the response of the 
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Table 1. Surface critical exponents and their definitions. 



surface magnetisation to a surface field, xii = ~^w^^ is given by Xii = xi^ = !)• From 
tlie profiles m{z) and one also obtains the surface excess quantities 

("^(^) - ^b) (5) 

2 = 1 

E (X(^) - X.) (6) 

2 = 1 

where m;, and Xb are the bulk magnetisation and the bulk susceptibility. Of further 
importance is the surface layer susceptibility, i.e. the response of the local surface 
magnetisation to a bulk field, usually denoted by xi- The critical surface pair correlation 
function behaves as 

where xi is the scaling dimension of the surface order parameter. For a (i-dimensional 
model p is a {d — l)-dimensional vector parallel to the surface. This behaviour of the 
surface correlation function differs from that of the bulk correlation function: 

Gfe(f-f') ~ |r-f'r'"' (8) 

where the value of the bulk scaling dimension Xb is usually smaller than that of Xi. Here 
r is a (i-dimensional vector. Critical correlations between a surface point and a bulk 
point exhibit a power law behaviour with the exponent xi + Xb- The surface critical 
exponents associated with the different critical quantities are listed in Table H 

The phase diagram of the three-dimensional semi-infinite Ising model is well 
established. The global phase diagram of the semi-infinite system depends not only 
on the values of the coupling constants and on the temperature but also on the external 
fields H and Hi ^UJ- I discuss in the following only the case of vanishing external 
fields. The qualitative results obtained in mean-field approximation [TTl [T21 ESI IHl U] 
represent the correct phase diagram as derived from renormalization group calculations 
[m Unj and from Monte Carlo simulations [131 HE] quite well, see Figure El If the ratio 
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Figure 2. Surface phase diagram of the semi-infinite three-dimensional Ising model. 



of the surface coupling Jg to the bulk coupling J^, r = Js/Jb, is sufficiently small, the 
system undergoes at the bulk critical temperature an ordinary transition, with the 
bulk and surface ordering occurring at the same temperature. Beyond a critical ratio, 
r > r^p ~ 1.50 for the semi-infinite Ising model on the simple cubic lattice [THl IT71 ITH], 
the surface orders at the so-called surface transition at a temperature T, > Tc, followed 
by the extraordinary transition of the bulk at T^. At the critical ratio r^p, one encounters 
the multicritical special transition point, with critical surface properties deviating from 
those at the ordinary transition and those at the surface transition. In mean-field 
approximation the critical ratio at the special transition point is = 1.25 for the 
simple cubic lattice. 

In the field theoretical formulation of the problem, see Eq. (0]), three stable 
renormalisation-group fixed-point values for cq can be shown to exist in absence of 
external fields. For cq = +oo resp. cq = — oo one is dealing with the critical behaviour of 
the ordinary resp. extraordinary transition. The fixed point corresponding to the special 
transition is located at cq = csp. Usually, a new variable c ~ cq — csp is introduced, 
such that the ordinary transition occurs for c > 0, whereas the extraordinary transition 
is observed for supercritical enhancements c < 0, the special transition point then being 
located at c = 0. 

Mean-field approximation yields similar surface phase diagrams for all 0{n) models 
in all dimensions. Of course, this scenario with three lines of continuous phase transitions 
is only realized in systems where the surface can sustain long-range order, independently 
of the bulk. Rigorous results ^Hj show that only the ordinary transition is encountered 
in the two-dimensional Ising model with short-range interactions in absence of external 
fields. This is readily understood as the surface is then one-dimensional. For continuous 
spins with n > 3 and isotropic couplings, a surface transition, where only the surface 
orders, does not exist in three dimensions, in accordance with the well known Mermin- 
Wagner theorem [SO]- Indeed, the surface effectively decouples from the bulk for very 
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strong surface couplings Jg ^ Jb and then forms an isolated two-dimensional system. 
An interesting case is that of the XY {n = 2) model which has a Kosterlitz-Thouless 
transition in two dimensions |2I1 122] • A phase diagram similar to Figure El has been 
shown to exist for the three-dimensional semi-infinite XY model [2S1 1^ 125! with the 
surface undergoing a Kosterlitz-Thouless transition for coupling ratios larger then some 
critical value, the bulk remaining disordered. Due to the peculiar properties of this 
transition, the multicritical special transition point is then of a special nature. Diehl 
and Eisenriegler have shown that in semi-infinite 0{n) models (and therefore also in the 
experimental relevant case of the semi-infinite Heisenberg model in three dimensions) 
anisotropic surface couplings may lead to an anisotropic special transition point and to 
a surface transition where the bulk remains disordered (2^1 ■ Phase diagrams similar to 
Figure 121 are found for all 0{n) models with isotropic couplings above the upper critical 
dimension. 

The models mentioned so far only cover a fraction of the research on surface critical 
phenomena. Other semi-infinite models studied include, for example, spin-1 Ising models 
|27j . Blume-Emery-Griffiths models j2Hl; Potts models ferrimagnetic Ising models 
IHO] , layered magnetic systems or bond percolation in the semi- infinite system (which 
can be regarded as the (g — > l)-limit of the g-state Potts model) IHSl ^^ 1^ IHUj. 

Before discussing the various surface universality classes, I want to mention two 
further interesting cases. For systems with a discontinuous phase transition, the surface 
order parameter may change continuously as the bulk transition point is approached. 
As discussed by Lipowsky |SZ] (see also |HH1) universal surface properties may show up 
for this case too. In the presence of this so-called surface-induced disorder, correlation 
lengths both parallel and perpendicular to the surface diverge at the bulk first-order 
transition point, thus inducing anisotropic power law behaviour for some bulk quantities 
|39j . A possible realization of this scenario may be encountered in the antiferromagnet 
UO2 where the surface layers order continuously while the bulk displays a discontinuous 
ordering |1U]. A second intriguing scenario is encountered in the g-state Potts model 
|4H 1321 where a phase transition to a low temperature phase with an ordered bulk but 
a disordered surface is observed. 

2.2. The surface universality classes 

It is obvious from the preceding discussion of possible surface phase diagrams that for 
every bulk universality class different surface universality classes may be realized. The 
different surface universality classes will be discussed in the following in more detail and 
their differences emphasized. 

At the ordinary transition both bulk and surface order at the bulk critical 
temperature. The ordinary transition is the only phase transition which is observed 
in the surface phase diagrams of all ferromagnetic three-dimensional 0{n) models in 
absence of symmetry-breaking fields. It has been studied intensively, both by analytical 
and by numerical methods. The surface critical exponents at the ordinary transition 
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can all be obtained by combining bulk exponents with one additional surface exponent 
Ai |121 113j . This new exponent results from the scaling function of the singular part of 
the surface free energy 

where t = {Tc — T)/Tc is the reduced temperature and hi is the surface scaling field 
which depends on both bulk and surface fields. The bulk exponent A;, is known from the 
singular part of the bulk free energy, whereas is the critical exponent of the excess 
specific heat Cg- Various scaling relations connect surface and bulk exponents. The 
critical exponents of excess quantities are obtained by combining bulk exponents: 



Ps = Pb- i^b (10) 

7s = 7b + (11) 

as = Ob - Vb (12) 
Other useful scaling relations are [d being the number of space dimensions) 

7s = 271 - 711 (13) 

7n = z/fe(t/-l-2xi) (14) 

71 = z/fc ((i - xi - Xfe) (15) 

(3i = vi,xi (16) 
Further scaling relations are discussed in [T31II1IS1- 



In Table E] estimates of various three-dimensional Ising surface critical exponents 
obtained by different techniques are compiled. One observes a very good agreement 
between the recent massive field-theoretical estimates and the numerical 

estimates. Similar agreement is also obtained for other 0{n) models. From the 
numerical point of view, the best investigated cases are n = ^Sl HEI and n = 1 
(Ising) im Ha CHI m El, whereas studies forn > 2 are scarce [21 UHl EOl EH- The 
rather few experimental determinations of surface critical exponents at the ordinary 
transition, using for example x-ray scattering at grazing angle, yield values which are 
found to be compatible with the theoretical estimates |52 j l58 t l5 H l55j. 

Exact results can by obtained by applying conformal invariance to semi-infinite 
systems. As shown by Cardy j2Z| conformal invariance yields in two-dimensional 
systems the exact values of the surface critical exponents at the ordinary transition and 
completely fixes the correlation functions. In higher dimensions it constrains the form of 
correlation functions near the free surface. Recent applications of conformal invariance 
include the determination of order parameter profiles in various two-dimensional systems 
with boundaries |SH] or the computation of three-dimensional Ising surface critical 
exponents from models defined on half sphero cylinders [HH] . 

The special transition point, located at the couphng ratio r^p, is a multicritical point 
where the bulk and the surface become critical. Three different critical lines (ordinary 
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Table 2. Estimates of Ising surface critical exponents in three dimensions at the 
ordinary transition, as obtained by different techniques. MF: mean-field theory, MC: 
Monte Carlo techniques, FT: field-theoretical methods. 









ill 


$ 


MF tlj 


1/2 


1 


1/2 


1/2 


MC gZI 


0.18(2) 


1.41(14) 


0.96(9) 


0.59(4) 


MC dHl 


0.237(5) 






0.461(15) 


MC 


0.2375(15) 


1.328(1) 


0.788(1) 




FT [60^ 


0.245 


1.43 


0.85 


0.68 


FT ^44! 


0.263 


1.302 


0.734 


0.539 



Table 3. Estimates of Ising surface critical exponents at the special transition point 
in three dimensions, as obtained by different techniques. MF: mean-field theory, MC: 
Monte Carlo techniques, FT: field-theoretical methods. 

transition, surface transition, extraordinary transition) merge at this point. The surface 
criticahty is thereby characterized by two new surface exponents: A^^ and $. The 
appropriate scahng ansatz for the singular part of the surface free energy then reads: 

fis^na) ^ |^|2-a. ^|^|-A, i^pAf |^|-* ^^^^ 

with the surface enhancement c ~ r — r^p. The crossover exponent $ governs the 
behaviour of the surface transition hne close to the special transition point: 

|T,(c)-T,|/T,~|c|^/*. (18) 

Estimates for various surface critical exponents obtained at the special transition point 
of the three-dimensional semi-infinite Ising model are gathered in Table 01 Again, good 
agreement between analytical and numerical estimates has been achieved. One also 
notes that the scaling laws (jl4|) - (jl6|) hold at this multicritical point. 

For surface enhancements exceeding the critical value Vgp two distinct phase 
transitions are encountered. The surface transition line, which is located at higher 
temperatures, separates the disordered high temperature phase from a phase where the 
surface alone is ordered. The bulk orders at the lower bulk critical temperature in the 
presence of an already ordered surface. This latter transition is coined extraordinary 
transition. The possible existence of a magnetic surface transition at temperatures 
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higher than the bulk critical temperature has been suggested to exist for various 
compounds as for example Gd [HIl El ESI El, Tb [HSl Ei, FeNig NiO [HHl Ei, 
NbSe2 i2Q| or Ni-Al alloys [H]. However, the experimental situation is usually not very 
clear. A good example for the encountered experimental difficulties is Gd. Whereas it 
was believed during many years that Gd undergoes a transition to a surface-ordered, 
bulk-disordered phase some 80 degrees above the bulk critical temperature, a recent 
study [72] claims that for pure Gd surface and bulk order at the same temperature, 
thus giving rise to the ordinary transition behaviour. 

Usually, the surface transition is considered to be of minor theoretical interest. This 
is due to the fact that at the surface transition the surface critical behaviour of a d- 
dimensional semi-infinite system is expected to be identical to that of the corresponding 
{d — l)-dimensional bulk system. Indeed, for perfect surfaces, the critical exponents of 
the two-dimensional bulk Ising model are encountered at the surface transition of the 
three-dimensional semi-infinite Ising model. However, the situation is not so simple for 
non-perfect surfaces presenting edges or extended surface defects. As will be discussed 
in detail in Sections 3 and 4, the local critical behaviour at non-perfect surfaces may 
be non-universal at the surface transition. In that case, local critical exponents which 
vary continuously as a function of the local couplings are encountered. Furthermore, the 
presence of the disordered bulk is refiected by the values of the local critical exponents. 
It is worth noting that these effects are not restricted to the Ising model but are also 
encountered in the three-dimensional easy-axis anisotropic Heisenberg model whose 
surface transition belongs to the universality class of the two-dimensional Ising model. 

At the extraordinary transition, the bulk orders in the presence of an ordered 
surface. The extraordinary transition with enhanced surface couplings and vanishing 
surface field is equivalent to the normal transition where the bulk orders in the presence 
of a field acting on the surface. This has been conjectured by Bray and Moore [73] and 
later proven by Burkhardt and Diehl for the Ising model [Z3]. The normal transition 
occurs in confined binary liquid mixtures at their bulk critical point. The equivalence 
of these two transitions is very useful for experimental studies as surface ordering 
fields are commonly encountered whereas physical systems with a genuine extraordinary 
transition are very scarce. To my knowledge, the extraordinary transition has only been 
investigated experimentally for NiO It should also be noted that some model 

systems possess a normal transition without exhibiting an extraordinary transition. 
A good example is the two-dimensional semi-infinite Ising model with short-range 
interactions. 

Surface critical exponents at the extraordinary transition can be completely 
expressed by bulk exponents. The singular part of the surface magnetisation, for 
example, varies close to the bulk critical temperature as m^'^"'^^ ~ |T — Tcf~°'''. 

To conclude this Section, let us finally mention the interesting possibility that 
the surface critical behaviour depends on the surface orientation. Examples where this 
dependence has been proven are binary alloys with a continuous phase transition [ZH] and 
Ising antiferromagnets in the presence of a magnetic field |75| . Ordinary transition 
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behaviour is encountered in these systems for symmetry-preserving surface orientations, 
whereas symmetry-breaking orientations lead to normal critical behaviour. 

2.3. Thm films 

Thin film magnetism is the subject of intensive current research activities, see [ZZj 
for a recent review. For the investigation of magnetic properties of thin films 
experimentalists have a large variety of experimental techniques at their disposal, 
ranging from ferromagnetic resonance to magneto-optic Kerr effect measurements. This 
has led to a large amount of interesting results concerning the magnetism of thin films. 

Critical phenomena in thin films have also been studied in recent years. In the 
following, I discuss some of the more general aspects of thin film critical behaviour, 
focusing thereby on thin films with perfect surfaces. The critical behaviour of more 
realistic films will be discussed in Section 4.2. 

First one has to remark that the critical temperature in thin films is a function 
of the number of layers forming the film. The temperature shift has been observed in 
numerous experimental studies and it has also been investigated extensively in theory 
[ZHl IZni IHOl IHH E2j- For thick films with L layers, the deviation from the bulk critical 
temperature Tc(cxo) is described by the well-known finite-size scaling relation [HH] 



The shift exponent A is given by A = l/z/?, where Ub is the correlation length critical 
exponent of the infinite system. In the ultrathin film limit, a linear dependence of 
Tc{L)/Tc{oo) on the film thickness is observed, see [HI] for a recent discussion. 

When analysing critical quantities in films with varying thicknesses, a dimensional 
crossover from three-dimensional criticality for thicker films to two-dimensional critical 
behaviour for ultrathin films is observed jHH IH3 EEl EI] , due to the truncation of the 
correlation length normal to the film |HH]. This crossover has been studied numerically 
in Ising films by analysing the critical exponent of the total film magnetisation |89j . 
Note that from a puristic point of view, the two-dimensional critical exponents should 
be observed for every finite film in a small temperature range near the critical point. 
However, as the width of this temperature window decreases rapidly for increasing film 
thickness, this temperature range may not be easily accessible in experiments or in 
computer simulations. Recent studies investigated this crossover in various systems 
using local effective critical exponents [^01 HH 112] • 

In thicker films of some systems, a crossover from a three-dimensional Heisenberg 
to a three-dimensional Ising behaviour may be observed jHZj- The formation of an easy- 
magnetisation axis is due to an increase of the magnetic anisotropy energy. In fact, 
the reduced symmetry at surfaces increases the anisotropy energy as compared to bulk 
systems where it usually is small. A change of the direction of the easy axis is often 
observed when changing the thickness of the film. For example, for Ni/Cu(001) the easy 
axis is in-plane for ultrathin films, but an out-of-plane easy axis is observed for films 
with more than 7 monolayers Long-range dipolar interactions contributing to the 




(19) 
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anisotropy are responsible for this behaviour. A further contribution to the magnetic 
anisotropy energy has its origin in the fact that (uhra)thin magnetic films are grown 
on a substrate. Indeed, a distortion of the lattice due to strain between the magnetic 
layers and the substrate may change the magnetic anisotropy energy as compared to 
the bulk system. The magnetic anisotropy energy is the subject of numerous theoretical 
studies, using different analytic methods [Hlj or ab initio techniques jHS]. In a statistical 
treatment of surface phenomena, the change in the magnetic anisotropic energy due 
to the presence of surfaces is usually modeled by effective short-range couplings with 
varying coupling constants, see [TH] for an example. Furthermore, when studying thin 
film critical behaviour theoretically, films are usually supposed to be free standing. 

Most thin films may be grouped into two different universality classes with respect 
to their critical behaviour. Films with an out-of-plane easy axis are theoretically 
described by Ising models (e.g. the Fe/Ag(100) system jHEj), whereas films with easy- 
plane magnetisation are modeled by XY-mode\s (e.g. the Ni/Cu(100) system ^Tj). 
Interestingly, a given material may exhibit both types of magnetisation, depending on 
the orientation of the film and the nature of the substrate. For example, Ni(OOl) grown 
on Cu(OOl) presents an out-of-plane easy axis, whereas for Ni(lll) on Re(OOOl) the 
magnetisation is in-plane. Spin reorientation transitions as a function of surface and 
bulk anisotropies have been studied theoretically for thin ferromagnetic films as well as 
for semi-infinite ferromagnetic systems in j98j . 

2.4- Surface critical behaviour near a Lifshitz point 

Competing interactions are encountered in a large variety of physical systems as, 
among others, magnetic systems, alloys or ferroelectrics ll()()| llOU 1102^ ll()3j . These 
interactions may lead to rich phase diagrams with a multitude of commensurately and 
incommensurately modulated phases as well as to special multicritical points called 
Lifshitz points. At a Lifshitz point, a disordered, a uniformly ordered and a periodically 
ordered phase become indistinguishable |lU4j . A large number of systems (as, e.g., 
magnets, ferroelectric liquid crystals, uniaxial ferroelectrics or block copolymers) have 
been shown to possess a Lifshitz point. 

From a theoretical point of view the best studied Lifshitz point is that encountered 
in the three-dimensional ANNNI model ll()()| llOlj . The Hamiltonian of this model, 
defined on a simple cubic lattice, reads 



with Sxyz = =tl- Here J > and k > are coupling constants. In the planes 
nearest neighbour spins are coupled ferromagnetically with the coupling constant J, 
whereas in z- or axial-direction competition between ferromagnetic nearest neighbour 
and antiferromagnetic next-nearest neighbour couplings takes place, leading to the 
appearance of spatially modulated phases. In the infinite system infinitely many 




xyz 



xyz 




(20) 
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commensurately and incommensurately modulated phases appear in the (T,k) phase 
diagram |106t I107t llUUj . In thin ANNNI films, however, only modulated phases 
compatible with the thickness of the film may be stabilized, leading to a different phase 
diagram for every film thickness |1U8| llU9t IllUj . Note that the related phase transitions 
in thin helimagnetic and incommensurately modulated films have also been the subject 
of recent studies |11H I112j 

The Lifshitz point encountered in the three-dimensional ANNNI model is of the 
uniaxial Ising type, the order parameter at this multicritical point having only one 
component. The term uniaxial denotes the fact that wave vector instabilities only show 
up in a single direction. In general, the Ginzburg-Landau free energy density in d 
dimensions may be written in the presence of a uniaxial Lifshitz point in the following 
form: 

/ = a2<^' + a4<^' + 6i |Vi + 62 

where is the order parameter, Vi the space derivative in the axial direction, whereas 
V(d-i) is the gradient operator in the directions perpendicular to that direction. The 
coefficient bi in equation ()21|) may change sign because of the competition between ferro- 
and antiferromagnetic interactions along the special direction. At the Lifshitz point, bi 
vanishes and the last term in equation 1)211) becomes relevant. 

A uniaxial Lifshitz point is only a special case of more general Lifshitz points |lU4j 
(see [11 3 j for a recent brief review) characterized by the number of space dimensions, d, 
the number of order parameter components, n, and the dimensionality m of the subspace 
where the wave vector instability occurs. The Lifshitz point encountered in the three- 
dimensional ANNNI model is then given by the set {d, n, m) = (3, 1, 1). Uniaxial Lifshitz 
points are strong anisotropic equilibrium critical points where the correlation lengths 
parallel and perpendicular to the axial axis diverge with different critical exponents 
and u^. The discovery of the existence of this kind of multicritical points in 1975 led to 
numerous theoretical studies of their bulk critical properties. With the exception of an 
early attempt by Gumbs |114j . surface critical phenomena at a bulk Lifshitz point have 
only been studied very recently |115| I116| 1117) 1118) 1119) ll2Uj . Most of the results have 
been obtained for the semi-infinite ANNNI model. 

Due to the anisotropy of the ANNNI model, surfaces with different orientations are 
not equivalent. The following two surface orientations have been considered (see Figure 
El): surfaces perpendicular to the axis of competing interactions (case A) and surfaces 
parallel to this axis (case B). As usual, the index b [s) indicates bulk (surface) couplings 
in the following. 

For case A modified surface couplings connecting neighbouring surface spins are 
introduced in addition to the usual ANNNI interactions, see Figure Et- Three different 
scenarios have to be distinguished, depending on the value of the bulk competing 
parameter Kb- When Kh is smaller than the Lifshitz point value = 0.27, the bulk 
undergoes a second order phase transition between the disordered high temperature 
phase and the ordered, ferromagnetic, low temperature phase at the critical temperature 



V(d_i) (j) 



+ Ci 
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Figure 3. Cross sections of semi-infinite three-dimensional ANNNI models showing 
two different types of surface orientations: (a) surfaces perpendicular to the axis 
of competing interactions, (b) surfaces parallel to this axis. Jt and Jg denote the 
nearest neighbour bulk and surface couplings, respectively, whereas the axial next- 
nearest neighbour interactions are labeled by the bulk, Kf,, and surface, Hg, competing 
parameters. Surface lattice sites are represented by filled points. 



Tc{K,b)- This phase transition belongs to the universahty class of the 3d Ising model. 
Consequently, the surface phase diagram will resemble that of the 3d semi-infinite Ising 
model. At the bulk Lifshitz point, Kb = i^bi recent mean-field treatment |115j 
yields for the semi-infinite ANNNI model a surface phase diagram similar to the Ising 
model, but with a set of critical exponents different from those of the Ising model. The 
predicted existence at the bulk Lifshitz point |115j of a special transition point is in 
variance with the earlier treatment of |114j but agrees with the Monte Carlo results of 
|118j . Finally, for axial next-nearest neighbour bulk couplings Kb > k^ the bulk phase 
transition from the disordered to the modulated phase belongs to the universality class 
of the 3D XY model |jl21j. It may then be argued that at the ordinary transition one 
is dealing with a critical semi-infinite 3D XY model with Ising-like surface exchange 
anisotropics. Surface exchange anisotropics being irrelevant near the ordinary transition 
PK] , one therefore expects the ordinary transition critical behaviour of the semi-infinite 
3D ANNNI model with Kb > k^ to be identical to that of the three-dimensional semi- 
infinite XY model [21] . A surface transition belonging to the 3D Ising universality class 
will again be encountered for strong surface enhancements. The ordinary transition and 
the surface transition lines are then expected to meet at an anisotropic special transition 
point |2E]- 

Case B consists of surfaces oriented parallel to the axial direction (see Figure Eb)- 
The introduction of modified nearest neighbour, Jg, and axial next-nearest neighbour 
couplings. Kg > 0, in the surface layer leads to intriguing and very complex situations. 
First, one has to note that the critical value of the coupling ratio r^p, needed for the 
surface to get critical by itself, depends both on the bulk, Kb, and on the surface, Ks, 
competing parameters. Three different cases may be distinguished, depending on the 
value of r = Jg/ Jb- 
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1) For r < Tsplnf,, k^) the ordinary surface transition is encountered, with an Ising 
bulk ordering (k^ < k^), a Lifshitz point bulk ordering = k,^), or a modulated 
bulk ordering (k^ > k^). 

2) For r > rsp{nh-, i^s) the surface orders at a higher temperature than the bulk. 
For Kg < 1/2 the surface transition belongs to the universality class of the two- 
dimensional Ising model, whereas for Kg > 1/2 a floating incommensurate phase 
appears in the surface layer. No surface transition, where the surface alone orders, 
will be encountered for Ks = 1/2, as in the 2d ANNNI model the paramagnetic 
phase extends down to T = for this value of |l(H)j . 

3) For r = rsp{nb, i^s) the special surface transition point is encountered where the 
ordinary transition line and the surface transition line merge. Depending on the 
values of Kf, and k^, very interesting possibilities arise. Whereas for < and 
Kg < 1/2 the usual scenario of an Ising ordinary transition meeting an Ising surface 
transition is encountered, for = k,^ and Kg < 1/2 a Lifshitz point ordinary 
transition merges with an Ising surface transition. These are the only cases studied 
so far in some detail |116| I118j . However, more exotic multicritical points may 
also be encountered. For example, one may have the case that at the special 
transition point an ordinary transition line with a modulated bulk (k^ > k^) meets 
a surface transition line to a floating incommensurate phase in the surface layer 
{Kg > 1/2). Finally, I should also mention the rather academic possibility that 
in four dimensions a Lifshitz point ordinary transition (/t5 = k^) merges with a 
Lifshitz point surface transition ( 

In their mean-field treatment Binder, Frisch, and Kimball |115| I116j considered 
surfaces oriented either perpendicular |115j or parallel |116j to the axis of competing 
interactions. For both cases the mean-field surface critical exponents at the Lifshitz 
point were determined at the ordinary transition and two different sets of critical 
exponents were obtained, see Tabled This dependence of the ordinary transition critical 
exponents on the surface orientation is explained by the fact that the bulk Lifshitz point 
is a strongly anisotropic equilibrium critical point. 

One should note that various scaling relations are fulfilled for case A |115j . as for 
example 



Here 7^, and are Lifshitz point bulk critical exponents which take in mean- field 
approximation the values (3^ = 1/2, 7^ = 1, and i^l^ = 1/4. Whereas (j2l|) also holds 
for case B, the scaling relations (j^^ and (j^^ have to be modified in this case. Indeed, 
the behaviour of excess quantities is governed close to a bulk critical point by the bulk 
correlation length along the direction perpendicular to the surface. The Lifshitz point 
being an anisotropic critical point characterized by two correlation lengths diverging 



Is = it + ^\\ 
Ps = Pt - 

7f = 27f - 7n 



(22) 
(23) 
(24) 
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/3f 


7f 




/5f 


7f 


case A 


OT/MF m] 


1 


1/2 


-1/4 


1/4 


5/4 




OT/MC [118 




0.62(1) 


0.84(5) 


-0.06(2) 


-0.14(4) 


1.69(7) 




SP/MC [118 




0.22(2) 


1.28(8) 


0.76(5) 






case B 


OT/MF 116 




1 


1/2 


-1/2 





3/2 




OT/MC [118 




0.687(5) 


0.82(4) 


-0.29(6) 


-0.46(3) 


1.98(8) 




OT/RNG [120J 


0.697 


0.947 


-0.212 


-0.462 


2.106 




SP/MC [118 




0.23(1) 


1.30(6) 


0.72(4) 







Table 4. Surface critical exponents at a uniaxial bulk Lifshitz point with the surface 
layer oriented perpendicular (case A) and parallel (case B) to the axial direction. OT: 
ordinary transition, SP: special transition point, MF: mean-field theory, MC: Monte 
Carlo simulations, RNG: renormalization group methods. 

with different critical exponents, should be used in case B instead of z/||". This then 
leads to the scaling relations 

f^s = f^t - (25) 

and 

Is = it + ^i- (26) 
Recent Monte Carlo simulations |118j revealed that the predictions of mean-field 
theory |115|I116] are qualitatively correct. As shown in Table 0] the values of the surface 
critical exponents at the ordinary transition in vicinity of the Lifshitz point indeed 
depend on the surface orientation. It is worth noting that for both surface orientations 
the value (3^ for the surface order parameter is clearly smaller than the corresponding 
value obtained in the semi-infinite Ising model. Mean-field theory yields for all these 
cases the same value [3f^^ = 1. 

Estimates for the Lifshitz point surface critical behaviour at the special transition 
point have also been included in Table lU Interestingly, the values of the critical 
exponents are very similar for the two different surface orientations. This is a strong 
indication that at the bulk Lifshitz point the surface critical behaviour at the special 
transition point may not depend on the orientation of the surface with respect to the 
axial axis. One may relate this observation to the fact that both the bulk (diverging bulk 
correlation length) and the surface (diverging correlation length for correlations along 
the surface layer) are critical at the special transition point. It may then be argued that 
the surface critical behaviour is governed to a large extend by the critical fluctuations 
along the surface, so that the surface orientation with respect to the direction of 
competing interactions is only of minor importance. At the ordinary transition, however, 
the surface is not critical and the surface critical behaviour is governed exclusively by the 
critical bulk fluctuations, leading to orientation dependent critical exponents because of 
the anisotropic scaling at the bulk Lifshitz point. 

Very recently Diehl and coworkers |119| I120j analysed the surface critical behaviour 
at bulk Lifshitz points using renormalization group methods. They thereby considered 
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general m-axial Lifshitz points where the wave vector instabihty takes place in an m- 
dimensional subspace of the ci-dimensional space. Thus for m = 1 one recovers the 
situation encountered in the ANNNI model. Restricting themselves to surfaces parallel 
to the modulation axes (i.e. to case B), they constructed the appropriate continuum 
101^ models and computed the critical exponents at the ordinary transition to order e^. 
Their results for m = 1 are included in Tabled One notes a very good agreement with 
the Monte Carlo results obtained in |118j . 

2.5. Critical dynamics at surfaces 

Besides changing the local static critical behaviour surfaces have also an effect on the 
dynamic critical behaviour [1221 [123 [123 1123 [121 [123 [123 [123 [133 [133 [133 [133 ITMl 
1135^ lllffif \1'A7\ llH8j . A central aspect of the works on dynamic surface critical behaviour 
concerns the possible classification of the distinct surface dynamic universality classes, 
similar to what has been done in the past for the dynamical bulk critical behaviour 
|lH9j . In that context semi-infinite extensions of the well-known bulk stochastic models 
are considered. Interestingly, different surface dynamic universality classes may be 
encountered for a given bulk model. This has especially been studied in the semi- 
infinite extension of the relaxation model B where the bulk order parameter is conserved 
113111132] . It has been shown that the presence of nonconservative surface terms, leading 
to a nonconserved local order parameter in the vicinity of the surface, yields a different 
dynamic critical universality class as compared to the case where these terms are absent. 
These two universality classes share the same critical exponents but are characterised 
by different scaling functions of dynamic surface susceptibilities. 

It is important to note that no genuine dynamic surface exponent exists |12l-{j . 
Indeed all exponents describing the equilibrium critical behaviour of dynamic quantities 
can be expressed entirely in terms of static bulk and surface exponents and the dynamic 
bulk exponent z. For instance the dynamic spin-spin autocorrelation function decays for 
long times as ^-^xi/z -^j^gj^g xi is the surface scaling dimension. The value of the dynamic 
exponent z is approximately 2.17 (2.04) in the two-dimensional (three-dimensional) Ising 
model. 

The effect of surfaces on nonequilibrium dynamics after a quench from high 
temperatures, T ^ T^, to the critical temperature has been investigated in \l'S'S\ I134| 
I138j . These studies are to some extent complementary to the investigations of universal 
short-time behaviour in the bulk |14()j . Indeed, similar to the bulk magnetisation, the 
surface magnetisation displays at early times a power-law behaviour 

mi(t)~mi,ot^ (27) 

with 6 = {xi — Xi)/z. Here mi^o is the small surface magnetisation of the initial state, 
Xi is the scaling dimension of the surface magnetisation, whereas the nonequilibrium 
exponent Xi is the scaling dimension of the initial magnetisation |14()j . A corresponding 
power-law behaviour is also obtained for the local nonequilibrium autocorrelation in 
the long time limit t ^ 1 |133^ I134j . It has been revealed recently |138j that in the 
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case Xi < Xi a new effect, called cluster dissolution, takes place, which leads to an 
unconventional, stretched exponential dependence of the short-time autocorrelation. 
A crossover to a power-law behaviour is then observed at later times. Interestingly, 
this stretched exponential behaviour is observed in the important case of the three- 
dimensional semi-infinite Ising model where Xi = 0.53 and xi = 1.26 |138j . 

There has been some progress in the analysis of dynamics in semi-infinite systems, 
but the situation is still very unsatisfactory as the effects of surfaces on dynamics have 
up to now only been studied in a very unsystematic way. From the experimental point 
of view the situation is even worse as no experimental studies of dynamic surface critical 
behaviour have yet been published. It has to be noted that some possible experiments, 
involving inelastic magnetic scattering of neutrons |141j or Mossbauer and NMR spectra 
|142j . have been discussed in the literature, but they have not yet been realised. 

2. 6. Surfaces and fields 

Up to this point we have only discussed surface phase diagrams in absence of external 
fields. In fact, the global phase diagram not only depends on temperature and on 
surface enhancement, but also on symmetry-breaking bulk and surface fields |10j. A 
discussion of the rich field of wetting phenomena, including for example critical wetting 
or prewetting, is beyond the scope of this article, even so some recent investigations of 
wetting criticality in wedge-shaped geometries will briefly be mentioned in Section 3.3. 
The reader interested in this field is referred to the review by Dietrich |143j and to the 
overview of experiments by Bonn and Ross |144j . I also refrain from dealing with the 
related topic of localisation-delocalisation transitions observed for example in Ising films 
with competing surface fields. This topic was the subject of a recent review article by 
Binder et al. 

Instead, I focus here on the normal transition and, especially, on the crossover 
between ordinary and normal transitions in presence of weak surface fields. This 
crossover has attracted much interest in the past and is also of relevance when discussing 
recent experimental studies of surface critical behaviour. 

At the normal transition the bulk orders at the bulk critical temperature in presence 
of an ordered surface. The non-zero surface magnetisation is thereby generated by a 
symmetry-breaking surface field. A similar situation is encountered at the extraordinary 
transition where the finite surface magnetisation at the bulk critical temperature is due 
to strong enhancement of the surface couplings. This two different transitions have been 
shown mj to belong to the same surface universality class. 

It seems natural to expect in the vicinity of the normal transition a monotonic 
decaying magnetisation profile for T > T^. This is indeed observed for strong surface 
fields, but for weak fields a more complex behaviour with a non-monotonic profile 
may be observed. Assuming that the surface enhancement c is subcritical, i.e. that 
without external fields one would be at the ordinary transition, it has been shown 
|146U147lll48j that a small surface field leads to a short-distance increase of m{z). This 
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critical adsorption in systems with weak surface fields has been studied subsequently in 
a series of papers |149| \lbl\ 11521 1153j . The increase of the layer magnetisation is 
due to the fact that a weak surface field gives rise to a macroscopic length scale, yielding 
a power law behaviour 



where Xb and xi are the bulk and surface scaling dimensions, respectively. This increase 
continues up to a distance 1°'''^ ~ f^i/{d~i-xi) ^ where a crossover to the normal monotonic 
decrease m{z) ~ z~^'' sets in. Recent studies of critical adsorption in a weak surface 
field for a homologous series of critical liquid mixtures |154j have permitted to observe 
the predicted increase experimentally. It is worth noting that for surface enhancement 
corresponding to the special transition point m{z) displays in the weak surface field 
limit a monotonic decay, but with two different power laws in the limits z — >■ and 
z — > oo |155j . 

It has been suggested |148| I156j that this crossover between the ordinary and the 
normal transitions may explain some puzzling experiments on surface critical behaviour. 
In these experiments |S^ll57j critical exponents compatible with the ordinary transition 
critical behaviour were measured, but at the same time the existence of residual long- 
range surface order was revealed at temperatures above the bulk critical temperature. 
This behaviour is readily explained by assuming that a weak surface field exists, yielding 
a surface structure factor governed by the ordinary behaviour in the case that 1°'^'^ exceeds 
the bulk correlation length. This interpretation assumes the existence of a surface field. 
It is therefore important to note that this kind of weak surface field can indeed arise in 
the studied compounds due to non-ideal stoichiometry effects |158[ I7H] . 

3. Surfaces with edges and corners 

The semi-infinite model with a fiat surface may be considered to be a special case 
of a more complex wedge geometry where two planes meeting at an angle 6 form 
an infinite edge. For 6 = it the fiat surface is recovered. Cardy showed that at the 
ordinary transition edge critical exponents depending continuously on 6 arise on purely 
geometrical grounds For a given opening angle 6, however, the values of the critical 

exponents are expected to be universal and independent of microscopic details such as 
the strengths of the coupling constants or the lattice type. Whereas edge singularities at 
the ordinary transition have been studied intensively, especially in two dimensions |161j . 
edge critical behaviour at the surface transition and at the normal transition have in 
general been overlooked. The recent investigations of edges and corners in systems with 
enhanced surface couplings and/or surface fields revealed some unexpected phenomena 
which will be reviewed in Sections 3.2 and 3.3. 
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3.1. Ordinary transition 

Cardy considered in his work d- dimensional 0{n) models containing edges formed by 
{d — 1) -dimensional hyperplanes meeting at an angle 9 [160J. For this geometry local 
critical exponents changing continuously with the angle 6 are already obtained in mean- 
field approximation. The edge magnetisation critical exponent for example is given 
by 

For the opening angle 6 = tt, the surface critical exponents are recovered. In the 
following, we use the same notation as Cardy and refer to edge quantities by the 
subscript 2. Corner quantities resulting, in dimensions d > 3, from the meeting of 
three hyperplanes will carry the subscript 3. 

Near the bulk critical point Tc, the edge energy density has the scaling form |16Uj 

f^ = \t\^^-^>''^[h\t\-^>',h\t\-^\h2\t\-^') (30) 

with t = (Tc — T)/Tc. Ufy is the critical exponent of the bulk correlation length, h, hi, 
and h2 are magnetic fields. /i2 only acts on edge spins, whereas h acts on all spins and 
hi on all surface spins. and Ai are the bulk and surface exponents appearing in the 
singular part of the free energy density of a semi-infinite system. The new exponent 
A2, which has been computed in first order of an e = d — 4 expansion in |160j . changes 
continuously with the wedge angle 9. All edge critical exponents may be expressed by 
A2 together with bulk and surface critical exponents. From the renormalization group 
point of view, the angle dependence has its origin in the invariance of the edge under 
rescaling. This makes the opening angle a marginal variable and may therefore lead to 
angle dependent local critical exponents. 

Inspired by these results, subsequent work mainly studied the infiuence of "edges" in 
two dimensions. In the following only some selected results obtained in two dimensions 
are discussed, the interested reader is referred to the excellent review of Igloi, Peschel 
and Turban |161j for a more complete account. The "edge" then reduces to one point 
and forms a corner in a two-dimensional system. Based on numerical and analytical 
calculations of isotropic two-dimensional Ising models |162j , the critical exponent of the 
local edge magnetisation was postulated to be 

ft = Ye 

for 2d Ising models. A similar equation was proposed for anisotropic lattices, with the 
opening angle replaced by an effective angle depending on the ratio of the different 
couplings. Using the conformal transformation w = z^l'^ , which transforms the semi- 
infinite system in the z-plane to a wedge with opening angle 9 in the w-plane, a simple 
relation between edge and surface critical exponents in two dimensions was derived 
[SZl ll^2j . Thus, for the local magnetisation one obtains 

/32 = ^ /3i (32) 
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which is in accordance with equation (|3ip for the Ising modeL Other work focused on the 
temperature behaviour of the local magnetisation in two-dimensional Ising models with 
various opening angles and lattice types I164| I165[ Up to now a complete 

solution has only been obtained for the square lattice Ising model with 6 = 7r/2 
|l(i7| 11681 Further studies also investigated the influence of edges in other two- 

dimensional systems. Examples are polymers |17m I17t| . regarded as 0{n) models in 
the limit n — > 0, or Potts models |172j . 

A qualitatively different critical behaviour is observed in systems with parabolic 
shapes |173t I161j . These systems have the remarkable property that they are 
asymptotically narrower than wedges. In this geometry one does not observe the usual 
power laws at criticality but stretched exponentials. Close to the critical point the tip 
magnetisation vanishes like exp(— a t-^) with a, 6 > and t = (T^ - T)/Tc. 

For small three-dimensional systems, edges and corners may be expected to play a 
dominant role, for instance, in nanostructured materials. Studies to reproduce average 
properties of such small clusters of atoms have been performed. These include several 
Monte Carlo investigations (see, e.g., jl74l I175| 1176] ). Critical phenomena at edges in 
three-dimensional models have, however, been addressed only rarely. Most investigations 
of edge criticality focused on Ising models, but some studies of percolation |177tll78j and 
of polymers |17m I179j in three-dimensional wedges have also been published. Whereas 
most of these studies of edge critical properties investigated the case that at the bulk 
critical point both surfaces forming the edge undergo an ordinary transition, Wang 
et al. |18()j also computed critical exponents for the cases that at least one of the 
surfaces undergoes a special transition. In this renormalization group study, the authors 
limited themselves to the opening angle 6' = 7r/2. They also studied the corner critical 
behaviour of cubic systems, each surface having von Neumann (special transition) or 
Dirichlet (ordinary transition) boundary conditions. Saxena |181j discussed, based on 
renormalization group calculations, the possible phase diagrams of Ising models with 
an edge formed by two perpendicular surfaces. Similar methods were used by Larsson 
|182j in his study of Ising edge critical behaviour, limiting himself mainly to the ordinary 
transition. The edge exponent A2(6'), see equation (jHUjl . was computed for some angles 9 
by Guttmann and Torrie |17nj using high temperature series expansions. Based on these 
results, they proposed an expression for /S.2{9) supposed to be valid for all 6. Early Monte 
Carlo simulations of Ising models were done by Mon and co-workers in order to compute 
the edge exponent /S.2{6 = n/2) \lKi\ as well as the corresponding corner quantity 
A3 for a cube |184j . In these studies, only rather small systems were investigated. 
The critical free energies of Ising models with edges and corners were also studied 
|185| I186j . following a similar investigation in two-dimensional critical systems with 
corners |187j . 0{n) models with edges and corners were studied in the limit n — > 00 
in jl88j . The temperature dependence of the edge and corner magnetisations as well 
as related quantities were only investigated recently for three-dimensional Ising models 
in the case of equal surface and bulk couphngs |189[ using modern simulation 

methods. 
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Figure 4. Geometry of a model with (100) and (DIG) surfaces, i.e. edges with opening 
angles 9 = 7r/2. Jf, and Js are the bulk and surface couplings, respectively. Jg is the 
coupling between neighbouring edge spins, J^s the coupling between an edge spin and 
its neighbouring surface spin. 
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Figure 5. Geometry of a model with edges with opening angles 9 = tt/A and 9 = Stt/A. 
Shown is a cut through the crystal perpendicular to the edge direction. The full 
lines show the couplings between neighbouring spins, whereas the (110) surfaces are 
indicated by the dashed lines. 



To introduce edges in Ising magnets, periodic boundary conditions along one axis, 
the z-axis, are applied. The remaining four free surfaces of the crystal may be oriented 
in various ways leading to different opening angles 9 at the edges. As shown in Figure 
m pairs of (100) and (010) surfaces lead to four equivalent edges with opening angles 
6 = 7t/2. The intersections of (100) and (110) surfaces form two pairs of edges with 
9 = n/4: and 9 = 37r/4. This is illustrated in Figure El Besides the bulk coupling 
Jb and the surface coupling Js, further couplings may be introduced |189| I191j : two 
neighbouring edge spins interact with the edge coupling J^, whereas an edge spin is 
coupled to its neighbouring surface spin by the edge-surface interaction Jgs- When 
studying edge behaviour near the ordinary transition all couplings can be considered to 
have the same strength: Je = Jes = Js = Jb- 
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7r/4 


7r/2 


37r/4 


TT 


MF [innj 


2.50 


1.50 


1.17 


1.00 


RNG TO 


2.48 


1.39 


1.02 


0.84 


HTS [170^ 


2.30 


1.31 


0.98 


0.81 


MC |183i 




1.38(6) 






Mc mni 


2.3(1) 


1.28(4) 


0.96(2) 


0.80(1) 



Table 5. Predictions for the Ising edge magnetisation critical exponents from various 
methods for four different opening angles of the three-dimensional wedge. MF: mean- 
field approximation, RNG: renormalization group theory, HTS: high temperature series 
expansions, MC: Monte Carlo simulations. 

The values of the local critical exponent /?2 of the edge magnetisation obtained 
in various studies of Ising models are compiled in Table El The expected angle 
dependence of the local critical exponents is thereby nicely illustrated. The values 
given by mean-field theory are systematically too high, according to the predictions of 
the renormalization group |160j . These predictions, in turn, seem to be systematically 
too large, as suggested both by high temperature series expansions |17()j and by Monte 
Carlo simulations |189j . The last two methods yield results which are in close agreement 
with each other. The value of P2 for 6 = 7^/2 derived from the Monte Carlo study of 
Mon and Valles |183j differs significantly from the value obtained from the recent studies 
using cluster algorithm |189j , presumably due to the small systems investigated in that 
early study. 

In order to study the universality of the computed critical exponents, one may 
investigate, at a fixed opening angle, wedges which are formed by different pairs of 
surfaces. For example, comparing 9 = 7i/2 edges formed by (100) and (010) surfaces or 
(110) and (110) surfaces show that, on the one hand, the local magnetisations differ, but 
that, on the other hand, the values of the local critical exponents are invariant against 
rotation of the crystal |189j . Note that the invariance of boundary critical exponents 
against rotation of the crystal has been discussed and partly even proven for edges 
in two-dimensional Ising models |163| I161j . Similarly, modifying the strength of the 
couplings does not alter the values of the critical exponents at the ordinary transition 
|189j . Furthermore, changing the lattice type is not supposed to change the critical 
behaviour either. The independence of the values of the edge critical exponents on the 
lattice type has been studied in two dimensions |161j . 

Investigations of critical phenomena near corners in three-dimensional systems have 
up to now been limited to the special case where the three surfaces forming the corner are 
mutually perpendicular. In Monte Carlo simulations, Ising cubes on simple cubic lattices 
with free surfaces are considered |184| I189[ I192j . The corner magnetisation critical 
exponent is thereby determined to be /^a = 1.77(5) |192j . This value is significantly 
lower than the mean-field value = 2 |161j . Analytical results beyond mean-field 
are not yet available. 
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Recent simulations of water in hydrophobic pores also illustrated the effect of curved 
surfaces on the values of local critical exponents |193j . In that study Brovchenko et al. 
observed a sharp crossover for cylindrical surfaces, yielding an asymptotic value for the 
local order parameter larger than the value j3i ~ 0.82 measured for slitlike pores. The 
critical behaviour characteristic of the ordinary transition is expected to be observed 
in this case not too close to as one is dealing with the liquid-vapour transition near 
a weakly attractive surface, see the discussion in Section 2.6. Of course, the geometry 
is slightly different to the wedge-shaped geometries discussed so far in the literature, 
but nevertheless the kind of systems discussed in jl93j seems to be a good candidate to 
study the effect of (generalised) edges on the local critical behaviour experimentally in 
the future. 

3.2. Surface transition 

Edge and corner critical properties at the surface transition have been usually 
overlooked. This problem was briefly addressed in |182j and |181j . but no systematic 
study has been done until recently |19H I192j . As shown in the following, intriguing 
phenomena are observed at edges and corners when the surface orders whereas the bulk 
remains disordered. Especially, non-universal local critical exponents are encountered 
at the surface transition of crystals with edges and corners. As discussed in the Section 
3.1, edge and corner critical exponents at the ordinary transition depend on the opening 
angle but do not depend on microscopic details of the model as for example the values of 
the interactions, the lattice type or the orientation of the surfaces. The edge and corner 
behaviour at the surface transition is in marked contrast to this, as, for a fixed opening 
angle, local critical exponents change continuously with the strengths of the different 
couplings |19H 1192] . As we shall see, the values of edge and corner critical exponents 
at the surface transition also reflect the existence of the disordered bulk. 

At the surface transition of the three-dimensional semi-infinite Ising model with 
a perfect surface, the critical fluctuations are essentially two-dimensional. The surface 
critical exponents reflect this reduced dimensionality, e.g. /3i = /32d = 1/8. On the other 
hand, the edge presents a local perturbation, acting presumably like a line defect in a 
two-dimensional system. Simple one- dimensional defects in the two-dimensional Ising 
model have been studied intensively [Ml 1113 IMl IUHl 1113 M IMl M- It 
was shown that in the vicinity of these defects the local magnetic critical exponents are 
non- universal |195j . As these exact results provide a useful framework for discussing the 
numerical findings of |191j and |192j . I briefly review in the following the main results 
concerning the local critical behaviour near defect lines in two-dimensional Ising models. 

The plane Ising model with a defect line |194j is an interesting system displaying 
non-universal magnetic critical exponents. A simple analysis on the relevance 
of perturbations shows that in two-dimensional Ising models an energy-like one- 
dimensional perturbation is marginal, yielding continuously varying local critical 
exponents. This is a consequence of the fact that in the unperturbed system the 
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correlation length critical exponent u;, equals 1 and the scaling dimension of the surface 
magnetisation operator xi is equal to 1/2 |161j . The non- universal behaviour was 
proven for the first time by Bariev |195j . He analysed two types of defect lines: chain 
defects, where a column of perturbed couplings with strength Jch is considered, and 
ladder defects, where modified couplings of strength Ji connect spins belonging to two 
neighbouring columns. Bariev's exact results demonstrate the dependence of the local 
magnetisation critical exponent on the values of the defect coupling. For the ladder 
defect, the local critical exponent is 



with 



2 

Pi = — arctan^ {''^Tl (33) 



= tanh(— ^)/tanh(— ^), (34) 



whereas for the chain defect one obtains 



2 

Pch = arctan {kch) (35) 



= exp ( -2 ( - — ^ ) ) . (36) 



with 

J is the strength of the unperturbed interactions, whereas T2d is the critical temperature 
of the two-dimensional Ising model. For both cases, enhanced (reduced) defect couplings 
yield a lower (higher) local critical exponent as compared to the perfect two-dimensional 
Ising model, /?2d = 1/8. 

In the following, I discuss the influence of edges with opening angle 9 = 7r/2, see 
Figure HI on the local critical behaviour at the surface transition. Because edges are 
one-dimensional, and all couplings in the models are of short range, edge quantities 
only become singular at the temperature Tg where the surface orders. Near the surface 
transition, where the critical fluctuations are of two-dimensional character, the edge 
then acts like a defect line in an essentially two-dimensional bulk Ising model. The 
edge coupling Je corresponds to a chain-like defect, the edge-surface coupling Jgg to a 
ladder- type defect. The change in the topology at the edge compared to the surface 
amounts to a complicated ladder- type defect. 

The value of the critical exponent is non- universal, it varies continuously as a 
function of the coupling strength J^s as shown in Table IHl These findings may be related 
to the reported results for two-dimensional Ising models with a ladder defect. The ladder 
corresponds to the edge, and the ladder coupling Ji reflects not only the edge-surface 
interaction Jgs but also the reduced connectedness to bulk spins at the edge compared 
to the surfaces. For J^s = Je = Js and Jg = 2Jh, the critical exponent of the edge 
magnetisation has the value /?2 = 0.095(5), significantly lower than the critical exponent 
of the perfect two-dimensional Ising model. Comparing the local critical exponent f3i 
near a ladder defect, see equations (fH^ and (jHH), with P2, one may attribute an effective 
ladder coupling Ji^^ to the edge with J^^-^ > J {= Jg). This effective enhancement of 
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Jes J Si Je Js 


Jes Jsj'^'i Je Js 


Je Jsj'^'i Jes Js 




0.095(5) 


0.176(5) 


0.127(5) 


jeff 


1.22 J, 


0.74 J, 


0.99 Js 



Table 6. Ising edge critical exponents at the surface transition for systems with 
opening angles 9 = 7r/2 and Js = 2Ji,. J^^^ is the effective defect coupling of a ladder 
defect having the same critical exponent as the edge. 





Jes Js) Je Js 


Jes J sl'^1 Je Js 


Je Jsl'^1 Jes Js 




0.06(1) 


0.26(2) 


0.14(2) 


l^star 


0.082 


0.210 


0.128 



Table 7. Ising corner critical exponents /33 at the surface transition for systems with 
opening angles 9 = t: /2 and Js = 2Jf,. Pstar is the local critical exponent near a star 
defect where three ladder defects with defect couphngs Ji^^ , see Table |H| intersect. 

the couplings is due to the influence of the bulk spins. When lowering the coupling J^s 
while keeping and Jb constant, the value of /? increases, as expected from equation 
(1221), see Table IHl The weakening of the edge interaction Jg has a less pronounced 
impact on P2 than the weakening of Jes- The close agreement of P2 with in the case 
Je = Jes/'^ is fortuitous. It is due to a compensation of a reduction in following from 
the increase in the effective ladder coupling stemming from the edge topology, and of 
an enhancement in (32 following from the decrease in the strength of the edge coupling. 

At the surface transition of three-dimensional Ising models corner criticality 
deserves to be analysed as well. As edges are local perturbations acting similar to 
defect lines in two-dimensional models, the corners of a cube may be interpreted as 
intersection points of three defect lines. As shown in Table |7| one again obtains critical 
exponents changing continuously with the strengths of the different couplings. 

To explain these flndings, note that the corners are intersection points of edges 
and recall that at the surface transition the critical fluctuations are essentially two- 
dimensional. The critical exponent of the corner magnetisation, /^s, has been related 
to that of the magnetisation at the intersection of three semi-inflnite defect lines in 
the two-dimensional Ising model. The critical exponent (Sstar for star defects formed by 
three intersecting ladder defects has been calculated by Henkel et al. j203| I204j . One 
may assign to the edges an effective ladder coupling Jf-^-^ by setting (3i = (32- Inserting 
this effective coupling into the equation derived for a star defect formed by three ladder 
defects |208| I204j then gives (3 star- The comparison of [3^ with (3 star yields a satisfactory 
agreement, see Table Of course, a more reflned analysis would focus on the fact that 
the correct geometry of the present problem is not that of a plane but that of a cone 
with three defect lines meeting at an angle 7r/2 |2()5j (see |2()6j for a recent discussion 
of Ising models with conical singularities). Furthermore, the rather complicated nature 
of the edge as a simultaneously ladder- and chain-type defect line as well as the effect 
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of the bulk spins on the corner magnetisation should then also be taken into account. 
3. 3. Wedges and surface fields 

The critical behaviour near edges has also been studied recently at the normal transition 
|213j . The authors investigated critical adsorption near edges due to the presence of 
symmetry-breaking surface fields in a wedge. Besides studying the problem in mean- 
field approximation, they also presented some exact results for the two-dimensional case. 
Especially, the critical edge exponent j32 was discussed. It was shown that at the normal 
transition local critical exponents also depend continuously on the opening angle. 

Critical adsorption near edges is, however, only one of the many intriguing 
phenomena encountered in systems with both wedges and surface fields. In the case 
of adsorption of a fluid on a solid substrate with wedge geometry |214l I215j it is 
expected from thermodynamic arguments that the liquid fills the wedge completely 
at temperatures T > Tp, where the filling temperature Tp is lower than the wetting 
temperature T]y of a planar, but otherwise identical, wall |216j . In a series of recent 
papers Parry et al. \2n\ I218| I219| 1220^ I221j studied the corner filling transition in 
detail, focusing especially on fluctuation effects and on the universality classes of filling 
transitions. One of their prediction was that critical filling, i.e. the filling transition for 
T — > Tp , could be continuous even in cases where the related wetting transition is of 
first order. They also studied the divergence of various length scales associated with 
this phase transition and predicted, basing themselves on a fluctuation theory, that the 
interfacial height Iq from the bottom of the wedge should diverge as Iq ~ (Tp — T)~^/^ 

The wedge filling transition can be studied in Ising models with wedge geometries by 
applying surface fields which favor one of the two phases in the wedge j223| I224t I225t I226j . 
Exact results in the two-dimensional case |224j permit to establish the existence of this 
kind of transition. Recent numerical investigations |225t I226j revealed a very good 
agreement with the theoretical predictions of Parry et al. In addition, a new type 
of interface localisation-delocalisation transition was revealed in the three-dimensional 
double wedge forming a pore with a square cross section |226j . The critical exponents 
of this transition can be related to the critical exponents of the filling transition in a 
simple wedge. 

3.4- Non- equilibrium systems 

The influence of wedges on critical behaviour has not only been studied in equilibrium 
systems but also in non-equilibrium systems. Frojdh et al. \207\ introduced an edge 
into a directed percolation process and analysed the impact this edge has on the non- 
equilibrium phase transition observed in this system. This work generalized earlier 
studies of directed percolation in a semi-infinite geometry j2()8j to the wedge geometry. 
Angle dependent edge critical exponents were observed in this non-equilibrium case, too. 
Directed percolation processes have also been studied in two-dimensional parabolic-like 



Critical phenomena at perfect and non-perfect surfaces 



29 



systems with a free surface at y = ±:Cx^ |2(J9t l21Uj . For k <1/ z, z being the dynamical 
exponent, the surface shape is a relevant perturbation and the tip order parameter 
displays a stretched exponential behaviour. In the marginal case, k = 1/z, a non- 
universal local critical behaviour is again observed. 

Other non-equilibrium absorbing phase transitions have also been studied in semi- 
infinite systems recently |21HI212j . but these studies have not been extended to wedge- 
shaped geometries. 

4. Critical phenomena at non-perfect surfaces 

In the preceding Sections we discussed critical phenomena in systems with various 
geometries: semi-infinite systems, wedge-shaped systems, and cubes. All the results 
presented so far have in common that only idealized, perfect surfaces were considered. 
However, real surfaces are often naturally rough, as steps or islands occur during growth 
processes or result from the effect of erosion. Furthermore, methods from nanoscience 
permit to create artificial structures on top of films. Examples include lines of adatoms 
or regular arrangements of geometric structures. All these defects have some impact on 
magnetic surface quantities. 

In the following I review theoretical studies where surface critical phenomena in 
systems with various surface defects have been investigated. Section 4.1 deals with semi- 
infinite (i.e. bulk terminated) systems, whereas Section 4.2 is devoted to the infiuence 
of surface defects on the critical behaviour of (ultra)thin films. 

4.I. Semi-infinite systems with surface imperfections 

Geometric surface imperfections (e.g., islands or vacancies) and impurities may be 
stable on the time scale of magnetic phenomena and thus lead to quenched surface 
disorder |227j . Some experiments indicate an enhancement of disorder near surfaces, 
thus pointing to the possible realisation of quenched surface disorder in systems where 
bulk disorder is negligible. In a theoretical description, this kind of surface disorder is 
usually mimicked by random surface fields or by random surface couplings. 

Early studies mainly focused on the global phase diagram observed in systems with 
random surface couplings or random surface fields. In |228j dilute semi-infinite Ising 
models with bond and site dilution both in the bulk and at the surface were considered. 
Mean-field results obtained for semi-infinite Ising models with random surface and bulk 
fields were presented in |229j . The possible types of phase diagrams were established 
and the phase transitions as well as multicritical points were discussed. Phase diagrams 
obtained for semi-infinite transverse Ising models with random surface and bulk fields 
were analysed recently in |23Uj . 

In a series of papers, Kaneyoshi (see, e.g., \2'M\ \'2'A'2\ \2Xi\ ) studied semi- infinite 
Ising systems with an amorphous surface layer in detail. The amorphous surface can 
be mimicked by choosing randomly weak or strong nearest-neighbour ferromagnetic 
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couplings between surface spins. Sometimes, he also considered random couplings 
between surface spins and the underlying bulk spins. Using different analytical 
approaches, surface magnetic properties as, for example, the magnetisation of the surface 
layer, were analysed. The corresponding phase diagrams were derived from the results of 
an effective field theory which includes correlations. One interesting result obtained in 
these studies concerns the location of the special transition point. Kaneyoshi observed 
that the critical coupling ratio r^p is shifted to higher values in systems with a diluted 
surface and he conjectured that this is a common effect of surface amorphisation. This 
was confirmed by a Monte Carlo study |H] where the strength of surface bonds was chosen 
randomly between two different values Jsi and Js2, the ratio d = Jsi/Js2 measuring the 
degree of dilution. For d = 1/10 the special transition point is located at the mean 
coupling ratio Vgp = 1.7(0.1) for the simple cubic lattice, noticeably larger than the value 
obtained for the perfect surface r^p ^ 1.5, see Section 2.1. In order to understand this 
shift we recall that in the limit r = {Jsi + Js2)/i!^Jb) ^ 1 the surface effectively decouples 
from the underlying bulk and can be regarded as a two-dimensional system. However, 
it is well known that at a given mean coupling {J si + Js2)/2 the critical temperature of 
the two-dimensional Ising model is reduced by randomness j2Mj . Therefore, for fixed 
r, increasing randomness shifts the line of the surface transition to lower temperatures 
which in turn shifts the location of the special transition point to larger values of the 
coupling ratio. 

Sometimes, alloy surfaces are mimicked by semi-infinite surfaces with randomly 
decorated magnetic and nonmagnetic atoms located on the surface |235j . Again, the 
location of the special transition point is observed to shift as a function of the model 
parameters. Quantities having an impact on this shift are for example the concentration 
of the nonmagnetic surface atoms or the spin of the magnetic atoms. 

The investigations reviewed so far almost exclusively focused on global phase 
diagrams of semi-infinite systems with non-perfect surfaces. The different authors did 
not try to study in detail the influence of these imperfections on the surface critical 
behaviour. However, as a certain degree of imperfections is unavoidable when studying 
surfaces experimentally, it is very important to clarify whether the surface critical 
exponents are robust against surface imperfections. 

A first step in this direction was taken by Mon and Nightingale in their study of 
the influence of a random surface field on the surface critical behaviour of the semi- 
infinite Ising model |236j . This work was motivated by the determination of the surface 
order parameter critical exponent (3i in a study of wetting phenomena of binary mixtures 
consisting of a polar and a nonpolar liquid |237j . The value obtained for [3i in that study 
was compatible with the value encountered in the semi-infinite Ising model with a perfect 
surface, even so (5i described the vanishing of the surface order parameter in the presence 
of a random surface field with zero average strength. Using a Harris type criterion |238j . 
Mon and Nightingale conjectured that, at the ordinary transition, a random surface 
field is irrelevant for the surface critical behaviour of the Ising model. They verified 
their prediction with the aid of Monte Carlo simulations. Using renormalization group 
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techniques it was shown in |239j that surface bond-dilution is irrelevant for the semi- 
infinite Ising model, the universality classes of the different transitions being those of 
the pure system. 

In an extensive study, Diehl and Niisser |227j derived Harris type criteria for various 
types of quenched surface disorder with the aim to assess the relevance or irrelevance 
of these random imperfections on the pure system surface critical behaviour. In his 
original work, Harris |238j studied the stability of critical bulk systems in the presence 
of randomness. The well-known Harris criterion states that, in the weak-disorder limit, 
bond disorder is irrelevant for the bulk critical behaviour provided the specific heat 
critical exponent a is negative. Disorder is relevant for a > 0, whereas the case a = 
is marginal. Looking at the Ising model, one concludes, based on the Harris criterion, 
that in three dimensions disorder is relevant as a ~ 0.11, whereas in two dimensions one 
encounters the marginal case a = which has attracted much interest |24()j . One should 
note that the condition for irrelevance of disorder resulting from a Harris type criterion 
is a necessary but not a sufficient condition for the stability of the pure system's critical 
behaviour. Nevertheless, a Harris type criterion is usually a reliable indicator for the 
irrelevance of randomness. 

The Harris criterion generalized in |227j to surface critical behaviour states that 
short-range correlated disorder is relevant or irrelevant depending on whether some 
(surface) susceptibility, which depends on the kind of randomness under investigation, 
diverges or is finite at the critical point. For random surface fields the quantity of 
interest is the surface susceptibility xn of the pure system, whereas for random surface 
couplings the relevant susceptibility is the local specific heat Cu As these generalized 
susceptibilities have a singularity of the form 

Xn~|T-T,r^" (37) 

at the critical point, the criterion indicates that the disorder is relevant for positive Fn 
(i.e. 7ii or an), but that it is irrelevant for negative Fn. 

This criterion has been applied by Diehl and Niisser in various cases involving 
surfaces. Thus, the presence of random surface fields is expected to be irrelevant at 
the ordinary transition for all 0{n) models with n < 3 if the dimension > 2, in 
accordance with the results obtained by Mon and Nightingale |236j for the special case 
of the semi-infinite Ising model {n = 1). Interestingly, Feldman and Vinokur recently 
showed |241j that weak quenched surface disorder destroys bulk ordering in the case of a 
system with continuous symmetry (as the XY model), leading to a power law decay of 
correlation functions and therefore to the appearance of quasi-long-range order. In the 
two-dimensional semi-infinite Ising model the perturbation caused by the random surface 
field is marginal so that the criterion does not yield a definite answer. This case was 
studied subsequently in |242t I243| I244t 1245] where it was shown that the surface critical 
behaviour of the 2d Ising model is described by Ising critical exponents with logarithmic 
corrections to scaling. For the special transition point, the Harris like criterion of |227j 
indicates that random surface fields are relevant in dimensions d < 4. However, the 
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three-dimensional semi-infinite Ising model with quenched random surface fields does 
not present a special transition point nor an extraordinary transition, as the surface 
transition line, where the surface orders alone, does not exist anymore. This follows 
from the fact that at finite temperatures random bulk fields destroy long-range order 
in the two-dimensional bulk Ising system |246j . In the cases where the extraordinary 
transition still exists in the presence of random surface fields, the criterion predicts 
that the disorder is then irrelevant. For random surface couplings, the perturbation 
is found to be irrelevant both at the ordinary and at the extraordinary transitions. 
More interesting is the situation at the special transition point where early estimates of 
the critical exponent an in three dimensions yielded small negative jfOl 124 7j or small 
positive values [H], implying that short-range enhancement disorder is close to being 
relevant in three dimensions. Recently, new field-theoretical estimates also resulted in 
negative values for an ^HIIS]- The Monte Carlo results obtained for three-dimensional 
semi-infinite Ising models with random surface-bond disorder |H] are also compatibel 
with the irrelevance of random surface enhancement for the special transition point 
critical behaviour. However, as short-range random surface enhancement is close to 
being relevant, one might expect long-range correlated enhancement disorder to be 
relevant at the special transition point |227j . 

Diehl and Niisser also studied the impact of surface-enhancement disorder at the 
special transition of a bulk tricritical system. As in this case the Harris criterion did not 
permit a definite prediction, renormalization group techniques were used in a second 
paper |248j to clarify the situation. 

The robustness of the surface critical exponents at the ordinary transition against 
two types of surface imperfections was established in a Monte Carlo study |H1. The 
studied imperfections correspond to an amorphous surface, mimicked by random strong 
and weak couplings in the surface, and to a simple case of corrugation, due to the 
presence of a step of monoatomic height superimposed on a perfect surface. For 
the amorphous surface, the Monte Carlo results showed that the dilution, at a fixed 
temperature T < Tc, decreases the magnetisation at and near the surface. The effective 
critical exponent, derived from the magnetisation of the diluted surface, follows in 
the Ising case closely that of the perfect case, yielding the same asymptotic value 
(3i = 0.80(1). Randomness in the surface couplings is therefore irrelevant for the 
asymptotic behaviour of the surface magnetisation, and even of minor importance 
for the corrections to scaling 8J. Similarly, the estimate for surface susceptibility 
critical exponent 71 is compatibel with the one in the perfect case. These numerical 
findings provide support for the conjecture of |227j that surface enhancement disorder is 
irrelevant at the ordinary transition. The robustness of the critical exponent f3i against 
dilution observed in the numerical study was later proven by Diehl |249j who derived 
upper and lower bounds on the surface magnetisation of the diluted surface and showed 
in a rigorous way that /?i takes the same value in the diluted system as in the perfect 
system. 

In jS] the effect of corrugation on surface critical behaviour was studied in Ising 
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Figure 6. Geometry of a model with two straight steps of monoatomic height. Jb 
and Js are the ferromagnetic bulk and surface couplings, respectively. 

systems where an extra half layer in the form of a strip-like terrace of monoatomic 
height was superimposed on a perfect surface, as shown in Figure IHl For this geometry 
the magnetisation at the step-edge deviates most significantly from the magnetisation 
of the perfect surface. The local critical exponent describing the behaviour of the step- 
edge magnetisation was found to have the value 0.800(15), in agreement with that of the 
magnetisation of the flat surface. However, in contrast to the case of random couplings, 
corrections to scaling are here distinctively different from those of the perfect surface. 
Diehl also considered this type of imperfections in |249j and showed in a rigorous way 
that the critical exponent of the step-edge magnetisation is identical to that of the 
magnetisation of a perfect surface. Note that at the surface transition local critical 
magnetic exponents are expected to be non-universal close to the step-edge. Indeed, 
the step-edge should then act like a defect line in a system governed by two-dimensional 
critical fluctuations. This is supported by Monte Carlo simulations which yield for 
Js = 2Jh the value 0.33(2) for the step-edge magnetisation critical exponent ||9D| in the 
Ising model. 

The study of magnetic properties of rough, corrugated surfaces near criticality is 
a rather new and promising field. The case of one additional terrace on an otherwise 
perfect surface may be generalized to the more complex situation of vicinal surfaces 
with terraces separated by equally spaced monoatomic steps |250[ 12511 1252| I253j . In 
nature, rough surfaces often result from a growth process and strongly affect the surface 
magnetisation. Diffusion-limited growth results in a rough growth front following a 
Poisson distribution. This may be realized in simulations by considering a surface formed 
of columns of random heights |251j . see also |254j . Layer-by-layer growth, however, 
may result in films with a finite number of complete layers and one partially filled 
layer, thus yielding a different type of roughness |255j . Whereas the different kinds of 
roughnesses discussed so far have been shown to be (or are supposed to be) irrelevant 
for the surface critical behaviour at the ordinary transition, this is not always the case. 
In a recent interesting work, Hanke and Kardar |256t 125 7j showed that self-affine rough 
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surfaces may give rise to novel surface critical behaviour with surface critical exponents 
different from those of the perfect case. Similar results have previously been obtained 
in two-dimensional systems with fractal boundaries |258[ I259j yielding new multifractal 
boundary exponents. These results in two dimensions have been interpreted in terms of 
a scale-dependent distribution of opening angles of the fractal boundary |259j . The same 
analogy has been evoked in the discussion of the surface critical behaviour of self-affine 
surfaces in |256l 125 7j . 

Hanke and Kardar also studied critical correlations in the vicinity of a regularly 
patterned surface (as might for example result from a lithographic preparation of the 
surface). Using a perturbative calculation in the deformations in height from a flat 
surface they showed that the leading power law decay of the correlations is the same as 
for a flat surface, but with a modulated amplitude reflecting the shape of the surface 



Extended surface defects have been studied intensively in two dimensions. In the 
Hilhorst-van Leeuwen model |26()^ one considers a semi-infinite square lattice with 
inhomogenous couphngs (alternatively, the triangular lattice has also been investigated). 
In the direction parallel to the surface one has couplings with a constant strength Ji, 
whereas the strength of the couplings varies perpendicular to the surface as a function 
of the distance y to the surface: 



This extended perturbation is irrelevant for u > 1/vb = 1, yielding the same 
critical behaviour as the homogeneous semi- infinite system. For uj < 1, however, the 
perturbation is relevant and a spontaneous surface magnetisation is observed at the bulk 
critical point. The most interesting case is = 1, where the perturbation is marginal. 
Here a spontaneous surface magnetisation is again observed for values of A larger than 
some threshold A^, whereas for A < A^ the scaling dimension xi varies continuously as a 
function of A: Xi = |(1 — A/Ac). This intriguing behaviour has attracted much interest 
in the past |161j and has recently led to the discovery of an up to then unnoticed effect 
in the short-time non-equilibrium critical dynamics at surfaces |138j . 

Finally, before turning to thin films with surface imperfections, let us briefly 
mention some recent work concerning surface critical behaviour in an Ising model 
with quenched random defects in the bulk. In the three-dimensional bulk Ising model, 
quenched random defects are a relevant perturbation leading to modified critical bulk 
exponents |26HI262] . The surface critical behaviour of the three-dimensional Ising model 
with quenched bulk disorder has been investigated both at the ordinary |263t I264j 
and at the special transition point |2(i3^ EES], using different renormalization group 
techniques. The crossover behaviour between these two transitions was studied recently 
in |266j . Interestingly, modified surface critical behaviour is encountered independently 
whether quenched surface disorder is present or not. The corresponding marginal case 
d = 2 has also been studied in detail recently |267t 1268] , yielding Ising surface critical 
exponents with logarithmic corrections to scaling. Boundary critical behaviour of q- 
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state random Potts models hase been studied in |2ti9j for 3 < g < 8. It should be 
noted that the problem of bond percolation has also been investigated in semi-infinite 
geometries [S21 ESI EH ESI EEl • In this case distinct surface percolation transitions have 
been identified and the values of the local critical exponents have been computed. 

4-2. Thin films with surface imperfections 

In a series of papers Aarao Reis studied the dependence of physical quantities in two- 
dimensional Ising stripes |27Ul I27H I272[ I273j and three-dimensional Ising thin films 
|254| I255j on the surface roughness. Using transfer matrix methods, Reis investigated 
stripes where the column heights were chosen according to a Gaussian distribution 
with mean L {L: integer) and variance (AL)^/2 |27()^ I27t] . Two different cases were 
considered: AL constant and AL = L/Lq with some constant Lq. Whereas the 
roughness becomes unimportant for large L in the former case, in the latter case the 
roughness increases with the mean thickness. Reis paid special attention to the finite-size 
scaling and to the finite-size corrections due to randomness. In |272j he also considered 
the case of noninteger mean L and showed that the free energy displays an interesting 
oscillating behaviour as function of continuously changing L due to oscillations in the 
mean coordination number. Recently |273j . he completed his investigation of stripes 
of random width by studying roughness of the more general form AL ~ with 
< (3 < 1. He also considered correlated roughness by imposing the maximal height 
difference between neighbouring columns to be not larger than one lattice constant. 
Interestingly, these computations showed that the correlations had no systematic effect 
on the corrections to scaling. More relevant to the understanding of the critical 
behaviour of real films are Reis' studies of films with rough surfaces |254t I255j . Using 
Monte Carlo techniques, he studied the same kind of uncorrelated roughness as for the 
stripes: Gaussian distribution of thicknesses with integer mean L and AL constant for 
all L or AL = L/Lq |254j as well as noninteger mean L |255j . These studies were partly 
motivated by the observation that some film quantities depend on the growth condition, 
the different growth mechanism leading to different types of roughness. Reis studied 
the magnetic susceptibility, the specific heat, and the total magnetisation of films with 
different thicknesses. Especially, he showed that the simple equation (fnijl connecting 
the critical temperature of finite films with the bulk critical temperature does not hold 
for noninteger mean thicknesses L, as Tc{L) shows a convex behaviour between two 
consecutive integer values of L |255j . He also studied the critical behaviour of rough 
films and thereby observed that the considered types of disorder are irrelevant: as for 
flat films, the critical exponents of the two-dimensional Ising model are recovered in 
rough thin films. Up to now only rough films with uncorrelated roughness have been 
studied. It would be very interesting to investigate the influence of spatial correlated 
roughness on the critical thin film behaviour as well. 

The impact of additional regular structures, located on the surface, on the critical 
behaviour of Ising films is studied numerically in jHOl llU9j . The additional structures 
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consist of one or two adjacent lines formed by adatoms as well as straight steps of unit 
height. One may introduce local couplings with different strengths, similar to what has 
been done in Section 3.1 for the edges. For example, the strength of the bonds connecting 
two neighbouring adatoms in the defect line or the strength of the coupling between the 
adatoms and the underlying magnetic film may be varied. The main finding of is 
that non-universal critical behaviour is encountered in Ising films with additional surface 
defects: the local magnetic critical exponents close to the defects depend continuously 
on the local couplings as well as on the layer thickness. Interestingly, the presence 
of an additional line is already enough to change the local critical exponent of a one- 
layer system by a large amount. Indeed, in the case where all the couplings have equal 
strength one obtains the value Pi ~ 0.084 < 1/8 = j32d for the critical exponent of the 
local order parameter |90j . 

It is worthwhile mentioning that related experiments on films with additional 
adatom lines have been published |274[ I275| I276j . even if there the lines of adatoms 
were not connected magnetically to the film. The magnetism of lines of Ad adatoms on 
Ag surfaces has been the subject of some theoretical studies |277l 1278] . The magnetic 
behaviour of thin films with large terraces has also been analysed experimentally |279j . 

Increasing the film thickness further, the semi-infinite system with surface 
imperfections is reached at the end. The following two typical scenarios must be 
distinguished: (i) for Js/Jb < rgp, bulk and surface order at the same temperature 
(ordinary transition), whereas (ii) for Js/Jb > the surface orders at a higher 
temperature than the bulk (surface transition). Vgp is the critical coupling ratio of 
the multicritical point, see Section 2.1. Therefore, a completely different behaviour of 
local quantities near the surface imperfections is expected in both cases when varying 
the number of layers L. Numerical simulations of non-perfect films are in complete 
agreement with this expectation jHO]. For Js/Jb < ^sp the surface critical exponent 
f3i ~ 0.80 is obtained in the limit L — > oo everywhere on the surface, independently 
of any additional surface structure, as discussed already in Section 4.1. In a film 
the local magnetisation near the defect closely follows the local magnetisation of the 
corresponding semi-infinite system for low temperatures. At temperatures where the 
correlation length is comparable with the thickness of the film a crossover from a regime 
with isotropic three-dimensional fluctuations to a regime dominated by two-dimensional 
critical fluctuations takes place. Increasing the film thickness the crossover temperatures 
approaches the bulk critical temperature, yielding in the limit of the semi-infinite system, 
L — > oo, the value /3; = /5i ~ 0.80 for the critical exponent of the local magnetisation. 
Choosing Js/Jb > ^sp, non-universal local critical behaviour close to the adchain is 
expected at the surface transition even in the limit L — > oo. This situation is indeed 
comparable to that of an edge at the surface transition, the edge corresponding to an 
extended defect line as discussed in Section 3.2. The local magnetic critical exponents 
should therefore depend continuously on the local couplings at the additional line. This 
has been studied in jHO] by examining the case Js = 2Jb. The results obtained show that 
near an additional line f3i is affected both by the values of the local couplings and by the 
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Figure 7. Specific heat of an Ising film with equal couplings consisting of one layer 
plus half a layer. Two distinct peaks are observed corresponding to two distinct phase 
transitions taking place in this composite system. 



presence of additional bulk-like layers leading to non-universal local critical behaviour. 

Further results have been obtained for systems with a straight step on top of film. At 
the ordinary transition, a straight step on top of a semi-infinite system does not change 
the local critical exponents (8j. In thin films the critical behaviour is however quite 
different |nn |ll()9j . Introducing a straight step by adding half a layer of magnetic adatoms 
to the surface of the magnetic film one observes two sharp peaks in the specific heat 
PmilU9|f279j . see Figure[71 pointing to the existence of two different phase transitions. In 
fact, when considering a system with one half layer on top of a film one is dealing with a 
composite system displaying in the thermodynamic limit two distinct phase transitions 
at two different temperatures. Consider for simplicity the case of a single layer plus 
half a layer. One phase transition then takes place at the critical temperature of the 
2d Ising model, ksTclL = l)/Js = 2.269... and one at the critical temperature of the 
double layer, ksTc^L = 2)/Js = 3.207(3). The value of the critical exponent of the step 
magnetisation at the higher temperature phase transition is 1/2, i.e. is identical to the 
value of the surface critical exponent of the two-dimensional Ising model. At the lower 
critical temperature the single semi-infinite layer then orders in presence of ordered 
surface spins, the ordered double layer acting at that transition as a surface field. One 
is therefore dealing with the normal transition. The same scenario is expected to hold 
for finite films with an additional half layer with Tc(L + 1) > Tc(L), yielding at Tc(L + 1) 
the critical exponent 1/2 near the step edge, in accordance with the numerical findings 
of [HO] and |lU9j . Note that very recently |279j . two sharp susceptibilities peaks observed 
experimentally in thin films with large terraces have been interpreted as the signatures 
of two distinct phase transitions. 
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Thin films with an additional terrace covering half of the surface are only one 
example of composite systems displaying two distinct phase transitions. Further 
examples include two-dimensional Ising models with a defect column where the bonds 
differ on the two sides of the column j28Uj . Ising quantum chains where couplings and 
transverse fields differ in the two half chains |281j , or layered magnetic systems consisting 
of thin layers of coupled Ising spins with S = 1/2 and S = 1 |31j. In |282j . wetting 
phenomena were studied in models consisting of two semi-infinite systems with unequal 
critical temperatures connected by a defect plane. 

Finally, let us briefly mention that in the context of spin systems with continuous 
spin symmetry thin films with selected surface defects (amorphisation of the surface 
layer |283j or steps j284j ) have been considered in some cases, but no systematic study 
of the critical behaviour of this kind of systems has been done up to now. 

5. Concluding remarks 

Surface criticality has been the subject of intensive studies in the last thirty years. 
Thereby a large variety of different methods (analytical, numerical, and experimental) 
has been used, yielding a host of interesting results, as reviewed in this work. Critical 
phenomena at perfect surfaces are now in general well understood, at least when dealing 
with static critical quantities. This is not really the case for dynamic critical behaviour 
at surfaces for which a coherent picture has not yet emerged. A major problem in 
this context is the total lack of experimental studies on surface dynamic properties at 
criticality. The situation is also not very satisfactory from the theoretical point of view, 
as only selected results on equilibrium and non-equilibrium surface dynamic behaviour 
at criticality have been published. Clearly, this is one of the most important aspects of 
surface criticality that warrants more attention in the future. 

Wedge-shaped geometries, which can be viewed as generalisations of semi-infinite 
systems, have also been discussed in detail in this review. It is encouraging that the 
effect of curved surfaces on the local critical behaviour has been observed in simulations 
of liquid-vapour transitions near a weakly attractive surface. This may point to possible 
experimental systems where this kind of problems can be studied. There are indeed 
a vast number of theoretical predictions for this kind of geometry, and experimental 
investigations are therefore welcomed. 

In the last years the focus of research on surface criticality has somehow shifted, as 
the main emphasis has been on more realistic surfaces. The facts that real surfaces are 
usually rough, displaying a variety of different surface defects, and that experimental 
physicists can create artificial structures on top of a surface directly lead to the question 
whether these quantities have an impact on local critical behaviour. We have presented 
a comprehensive overview of the field, thereby discussing in detail critical behaviour in 
semi-infinite systems with surface defects as well as in thin films with additional surface 
structures. Some common surface defects have been shown to be irrelevant for the 
surface critical behaviour at the ordinary transition. There are however some interesting 
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exceptions, as for example the case of self-affine rough surfaces. The situation is even 
more complex at the surface transition where in a three-dimensional system the critical 
fluctuations are of two-dimensional nature. Indeed, additional structures as steps or 
lines of adatoms have been shown in numerical studies to lead to non-universal local 
critical behaviour where the values of the local critical exponents reflect the strengths 
of the coupling constants as well as the presence of the disordered bulk. It is worth 
noting that in thin films this kind of additional surface structures in general leads to 
non-universal critical behaviour. 

It is obvious from our overview that a large number of recent studies of the effects 
of surface defects on the local critical behaviour are either of purely numerical nature 
or are using rather crude approximations. This is especially the case when dealing 
with non-perfect surfaces at the surface transition. There is a need for more elaborate 
analytical approaches, and it is one of the intention of this review to encourage further 
theoretical (and also experimental) investigations of critical phenomena at non-perfect 
surfaces. 
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